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^ ■ The theory of recurrence relations of linear multi-component and multi- 

■ parameter systems on the basis of the canonical transformations theory of the 

dynamical systems' sets is constructed. The parameters of the grating's knots 
\ are defined from the condition of the invariance of the model under shifts along 

the grating. The connection with a zero curvature representation for models 

00 ' 

• on the grating is installed. The examples of two- and three-parameter sys- 

tems described by the hypergeometric functions M{a, [3, t) and M(a, (3, £, t) 
are considered in details. The canonical recurrence relations increasing and 
O ! decreasing parameters {a,/3,£} for solutions of the corresponding equations 

are constructed. 

43 I. INTRODUCTION 

> ■ 

a ■ 

Many equations of the theoretical physics especially of quantum mechanics are linear. In 
this connection there is a possibility of constructing solutions of these equations searching the 
canonical transformations of the linear dynamical systems (but this canonical transforma- 
tions are different from those of quantum theory). These transformations can be formulated 
in terms of the creation and annihilation operators for some "quantum numbers" of the 
corresponding quantities that conserve their values. The indicated possibility was consid- 
ered in H], where the method of constructing the recurrence relation for searching for the 
eigenvalues and eigenfunctions of linear operators is proposed. This method is restricted 
by the case of one- dimensional systems described by the Lagrangian containing the spectral 
parameters in a potential part only. 

In the present paper the generalization of this method for the case of the multidimensional 
systems and systems with the Lagrangian containing spectral parameters in a kinetic part 
too (for example, for the radial equation in hydrogen-like atom theory) is constructed. 

The physical meaning of the considered canonical transformations for the classical theory 
differs from that for the quantum theory. The general point is that these transformations 
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can be considered as translations on the multidimensional grating in the phase space. In 
the classical theory the parameters of the grating's knots define the different configurations 
of considered class of dynamical systems (for example, the set of mechanical oscillators with 
different masses and spring rigidities). Conditions of the invariance under the canonical 
transformations along the grating define the admissible grating parameters and, therefore, 
solutions and a spectrum of oscillations. 

In quantum theory the parameters of the knots give the possible states of the same 
quantum system. In its turn the translation along the grating is defined by the creation and 
annihilation operators stipulating the transition from one state to another. 

The canonical recurrence transformations obtained in this way are the analogies of the 
known Beclund's |§ transformations for finite-dimensional dynamical systems. Note that 
they do not coincide by their form with the standard recurrence relations cited in hand- 
books. The difference between them is explained by the way of their construction. The 
main here is that the canonical recurrence relations conserve the symplectic structure of the 
corresponding dynamical systems in contrast with the standard ones. That is why from 
the standpoint of the theory of dynamical systems the canonical recurrence relations have a 
certain advantage. This relations can be reduced to the standard ones by the corresponding 
normalization of eigenfunctions. 



II. A CANONICAL MAP 

Let us consider the sets V = {q,p,H} and V = {q,p,H} of the s-dimensional dy- 
namical systems with the generalized coordinates q = q(t) = .., q = q(t) = 
{q l {t) } . . . ,<f (t)}, momenta p = p{t) = {pi(t), . . . ,p a (t)}, p = p(t) = {pi(t), . . . ,p a (t)} and 
Hamiltonians H = H(p,q,t), H = H(p, q, t). 

We shall call the mapping T : T> — > T> the canonical one, if it is defined by the set 
T = {F} of generating functions F = F(q, q, t), so that for an arbitrary system (q,p, H) EV 
and for some function F G T there is the system (q,p, H) G V, for which the relations 

OF dF ~ dF 

Pi=Oj> Pl = ~W H = H + ^ (hJ,k,... = l,2,...,s), (2.1) 

are realized. 

If T> = V, we shall say that given set of dynamical systems T> is invariant under the 
canonical transformation F. 

Let D be a parametrized class of dynamical systems described by the Lagrangian 

L = L(q, q, t, A, fi a ) = - [g^q - (Uy - Xg^q 3 , (2.2) 

where q = dq/dt, g^ = gij(fi a ,t), Uij = Uij(fi a ,t) are some functions, det \ \gij\ \ ^ 0, {A,/i a } G 
D, D is some set of parameters, a = 1, . . . , d. The ordinary summation convention works 
here and on. 

If we write the Lagrange- Euler equations for the system ( |2.2| ) in the form 



Pk k = 9 H [~J t 9ikj t + q k = \q l 



(2.3) 
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we shall make a conclusion that A is a spectral parameter of the eigenvalue problem of the 
rf-parametric set of the linear operators Pl(fi a ,t). For the Hamiltonians of the set (|2.2f) we 
have 

1 .. l 

H(p, q, t, A, Ha) = -g tJ PiPj - -(Uij - Xgij)q l q\ (2.4) 

where 

Pi = 9ijQ j , (9 ik 9kj = (2.5) 

The condition of the invariance of the set T> of dynamical systems under the canonical 
transformations JF is 

H(p,q,t,X,fi a )=H(p,q,t,X,fi a ), W{X,fi a }eD. (2.6) 

Under sequential actions of the series of the canonical transformations a sequence of 
the systems of the type Q2.2| ) and the parameters {A m ,/i n Jm, n a = 1,2,...} G D appears. 
Introduce the collective indices K = (m, n a ), K = {m,h a ) and designate (A m ,/i na ) = ok-, 
(Am,/in a ) = &k- Let the corresponding coordinates and momenta be qK,qR an d Pk,Pk- 
Further, F K % = F K f c (qK,q k :t) be the generating function of the canonical transformation 
F: V -> V, so that {qK,PK, H(p K ,q K ,t,a K )} -> {q^,p k ,E(p k ,q k ,t,a k )). Then the 
conditions (pTTT) , (|2.6|) give the following equation for the generating function 



dF Kk + R ( ^ ^ ^ \ =H (_?FkK } qRi a . f t ) . (2J) 



dt \ dqx J \ dq 



K 



Further we shall write the indices K and K for the generating functions only, and for an 
arbitrary physical variable / we shall write / or / instead of fx or f k respectively. Then 
it is possible to rewrite the condition of the invariance ( |2.7|) for the system (|2.2|) with the 
Hamiltonian Q2.4|) in the following form: 



+ 9 3 dqi dqj - 9 3 -^T-q~T ~ ( U v - X ^)q Q J ~ (Uij - X 9ij )q q>. (2.8) 



III. A RECURRENCE RELATIONS AS A CANONICAL MAP 

Let us consider the problem of searching for the solutions of Eq . ( ^8|) ■ The canonical map, 
that converts the linear dynamical systems into linear ones, must be linear. That is why the 
generating function F Kk must be a quadratic function of the generalized coordinates q, q. 
Therefore we search for F Kk in the form 

Frk = \ (27^W - h 3 qy - atfq?) , (3.1) 

where 7^, by, are some functions. Substituting Eq.(|3.1|) for Eq.(j2.8|) and equating the 
coefficients at q l qi , q l q> , q l (p , we obtain 
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27rs - g l3 {lirbjs + IjAs) = ~§ K 1 {^isCjr + Ijskr), (3.2) 

-b rs + g tJ b ir b js + \g rs - U rs = g l3 i ir lj S , (3.3) 

- C rs + g l3 "1irlj s = g lJ C ir C js - U rs + A# rs . (3.4) 

If we find one of the particular solutions 7^ , bij, Qj, then, using Eqs. (|2.1| ), ( P-5|) , (|3.1| ), it will 
be easy to construct the binomial recurrence relations for the eigenfunctions q(t) and q(t) 
of the set of operators P l k 

(bij + gajA <? J = HjO 3 \ (3-5) 



z ii ~ 9ijj t ) Q J = lijQ 3 ■ (3-6) 



The obtained relations, increasing and decreasing K, are the analogies of the well-known 
Beclund's transformations || for finite-dimensional systems. They are also the generalization 
of the known relations of the factorization method 0, ||. It means, that Eqs.(|3.5|), ( p.6| ) 
install the correlation between the canonical transformation method and factorization one. 

The algebraic trinomial recurrence relations can be obtained considering the sum of the 
generating functions of the two sequential transformations 

F KkR = F KK + F k~K- ( 3 - 7 ) 

According to Eq. (|2.1| ) we have OF K ^^/ 'dq = 0, from which the linear algebraic recurrence 
relation 

lijQ 3 ~ (cij + hj)^ + y tj q = 0, (3.8) 

follows. Here 6^,5^,7^ are the coefficients of a quadratic form of the type (|3.1| ) for the 
generating function F^-.. 

The generating function F R ^ of the canonical transformation (q,p) — > {q,p), which 

is a composition of the canonical transformations (q,p) — > {q,p) and (q, p) — > {q,p), can 
be constructed excluding the intermediate state q from Eq.( |3.7| ) by means of Eq.(|3.8|). It 
makes possible to simplify the procedure of searching for the generating function of the 
transformation K — ► K decomposing it into more simple steps of the calculating the ele- 
mentary generating functions of the transformations: K = {m, m, . . .} — > {m+1, ni, ...}—> 
{m + + K, and after this to construct their composition. 

From the standpoint of classical mechanics we consider here the canonical transformation 
theory of an infinite-dimensional system with the action 



S = I dt HWfl l K^K + i U H ~ ^mgij)qKQK} 
K 



describing a multicomponent model on a (d + l)-dimensional grating. The values of the 
parameters ok = {A TO , // ni , . . . , fM nd } that define grating knots are determined from the 
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condition of an invariance of the dynamical system under grating translations induced by 
the canonical transformations of the system. 

One more interesting possibility connected with the method of the inverse scattering 
problem is shown below by the examples of more simple dynamical systems. 

For the particular case of the set of (i-component three-parametrical dynamical systems 
(s = 1) let us denote: gn = m{y, t), Un = m(i/,t)U(/j,,t), where fi, v are parameters. Write 
the Lagrangian (|2.2|) in the form 



m(u, t) 



q 2 + (U(^,t) - X)q 2 



It is convenient to find the generating function F K f> in such a form 



KK 



mS, S = —(2'yqq — bq 2 — cq 2 ) 
2a 



(3.9) 



(3.10) 



where 7 is some unknown constant, a, 6, c are desired functions of t. The conditions ( f2~T 
may be rewritten in the form 



dS ■ 

q = ~q j q 

aq 



m OS ~ TT d . n . 
H = H+—(mS). 



m dq 



The recurrence relations 



d 



b + a- q = ^q, 



' m d \ _ 
c- a— — g = 7g 
mat 



(3.11) 



(3.12) 



and equation for the function S 
2 d 



(mS) 



dS\ 2 m (dS\ 2 



-(x-uy--(x-u)q 2 

m \ at 1 m 



0. 



(3.13) 



mdt \dq J 

follow from Eq. (|2^) , ( |3.11| ). Eqs. (|3.2j )-(|3^) for determining a, b, c, 7 acquire the form 

mm 

a = — a H c — b, (3-14) 

m m 



-c + b= (X-X + U -U)a, 
m 



(3.15) 



771 

ab= (A-f/)a 2 + -(6c- 7 2 ). 
m 



(3.16) 



The system of equations ( |3.14|) -( j3.15|) is linear with respect to the unknown functions a, b, c 
and does not contain the constant 7, while the equation ( f3.16| ) is quadratic and contains 7. 
The procedure of the searching for the solutions of the system ( |3.14 )-( j3.16 ) consists in the 
finding a more simple particular solution of the homogeneous linear undetermined system 
( P-14| )- (|3.15| ) for which the quadratic equation (|3.16| ) is satisfied at the some unknown 
parameters {A,/i, 9}. Then the constant 7 will be determined. 
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The algebraic trinomial recurrence relations (without derivatives) and their compositions 
can be found by the way analogous to obtaining Eq. (|3.8|) (see ]l| also). As a result we shall 
have 

lKK+l (ck+IK+1 . b K +lK+l\ . 1K+IK+2 n /, <m 

Qk 1 qx+i H qx+2 = 0, (3.17) 

a,K \ &K CLK+l ) dK+l 

where K is the collective index of the totality of the parameters {A m ,//„, Ui}, i.e. if K = 
{m, n, 1} then K+ 1 means that one of the indices m, n, I is increased by one. The coefficients 
in Eq. (|3.17|) correspond to that of the generating functions Skk+i and Sk+ik+2 of the type 



( |3.10| ) of the transformations qx — ► qx+i and qx+i — * qx+2- 

Writing the similar equations for the sequence of the transformations qx+i — * qx+2 — * 
qx+3 and excluding the coordinate qx+2 by means of Eq.( |3.17| ), we obtain the two relations, 
which can be written in the matrix form 

Qk+2 = B K Q K , (3.18) 

where 

Qk=(" K ), Qk + 2=(1 K+2 ), (3.19) 



and Bk is a 2 x 2 matrix, which is expressed through the coefficients of the quadratic forms 
{Skk+i, ■ ■ ■ j Sk+2K+z\ °f the above sequence of the transformations; its explicit form is not 
indicated here because of its inconvenience. 

From the other hand the differential equation ( |3.12| ) can be rewritten in these terms as 

^ = A K Q K , (3.20) 
at 

where 

_ _ j ~ m K+li>K K m K+l1 K K+1 j 2^ 



The conditions of compatibility of Eqs.( p.l^ )- (|3.19| ) 

dB K 



dt 



+ B K A K - A K+2 B K = (3.22) 



are equations for the coefficients clk^ ^>k-, c-Ki Ik of the sequence of the generating functions 
Skk+i and represent zero curvature conditions for the models on the grating p|. They are 
the conditions of an applicability of the inverse scattering problem method for integrable 
models. In our case they are realized as a consequence of the equations of the type (|3.14 )- 



Thus the regular way of the constructing zero curvature conditions, which can be gener- 
alized for nonlinear dynamical integrable systems, follows from our analysis. 
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IV. THE CONFLUENT HYPERGEOMETRIC EQUATION 



Many equations of quantum mechanics under the separation of variables and after isola- 
tion of an angular part lead to the second order ordinary differential equations, which can be 
reduced to the hypergeometric (or confluent hypergeometric) equation. In this connection 
it is important to consider the canonical transformation theory of the equations of such a 
type in general case. 

Let us consider construction of the canonical recurrence relations for the confluent hy- 
pergeometric equation 

tq+((3-t)q-aq = 0, (4.1) 

where a, (3 are some parameters. It is the Lagrange-Euler equation for a dynamical system 
with the Lagrangian 

L = \t^-\e + - t d 2 )- (4.2) 
Comparing Eq. (|4.2|) with Eq. (|3.9|), we can conclude that 

A = 0, u = P, fi = a, U(a, t) = -, m((3, t) = lPe~ % . 

t 



Therefore the system ( [3.14j) -( 3.16|) has the form 



m m m ■ a — a 

a = — a H c — o, — c + b = a (4.3) 

m m m t 



a 2 m 



ab = a z + — (bc-Y) (4-4) 

t m 

According to the above-mentioned, we shall find the generating functions of the elementary 
canonical transformations a) {a,/?} — > {a, (3}, b) {a, /?} — >■ {a, 0} only. Besides note that 
Eqs.( |4.3| )-(0) can be written as equations with rational on t coefficients. Therefore it is 



naturally to find the solution in a class of rational functions. 

A. A canonical transformation {a,/?} — > 
In this case (3 = (3, m = m and the system ( [4.3| )-( [4~4"D acquires the form 



\ • a — a 
d = ( l)a + c — b, c + b= a, (4.5) 



ab = bc-rf - a 2 -. (4.6) 



t 

We shall obtain the simplest rational solution of the undetermined system ( |4.5| ) by assump- 
tion a = t. From ( f4.5|) we shall find 
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b = -(a — a — l)t + b , c= -(a - a + l)t + b Q - j3 + 1, 



(4.7) 



where 60 is the integration constant. Substituting Eq. (f4.7| ) for Eq. (|4.6|) and equating the 
coefficients at the same powers of t, one obtains 

b = a, a = a + 1, 7 = \Ja{a — (3 + 1) (4.8) 

We have restricted ourselves here by the positive value of the radical. As a result we have 

1 

'a,/3\a+l,l3 - 2' 



2Ja(a - (3 + l)qq - aq- (t + a - (3 + l)q 



-2 



(4.9) 



for the generating function ( |3.10| ). The canonical recurrence transformations (13.12) have the 
form 



d 



d 



a j r t—]q = ^a(a-f3 + l)q, [ t + a - (3 + 1 - t— ] q = y/a(a - (3 + l)q. (4.10) 



dt 



To determine the connection with the standard recurrence relations and thus with the con- 
fluent hypergeometric function M(a,(3, t), we shall make the following substitutions 

q = y{a)^N{^), q = y(a) ^N{a + 1), y(a) = M{a, (3, t) (4.11) 

in Eq.( 4.10p . From the condition of coincidence of the obtained recurrence relations with 
the standard ones |4[] (rather with their consequence) 



a + t— j y(a) = ay(a + 1) 



a - (3 + 1 - 1 - t— J y(a + 1) = (a - /3 + l)y{a) (4.12) 



the functional equation for the normalization multiplier N(a) 

aN(a) = (a-/3+ l)N(a + 1) 



(4.13) 



follows. It can be solved by using the gamma-function's property T(z + 1) = zT(z) from 
where the substitutions 



a 



r(a+l) r(a-/3 + 2) 

a — f3 + 1 — 



(4.14) 



(4.15) 



r(a) ' " r ' " r(a-/3 + l)' 
follow. Using these substitutions, equation ( |4.13|) can be reduced to the form 

Tja-p + l) A „ T(a-f3 + 2) r(q-/9 + 3) 

tv~^ ^ " = — rv — Zri — ( + 1 = fv 7 n ^ a + 2) = ■ ■ ■ . 

1 (a) 1 (« + 1) 1 (a + 2) 

The additional equations arise owing to considering the next elementary steps by the pa- 
rameter a. From here it is easy to see that T(a — (3 + l)N(a)/T(a) — G\ = const. Therefore 

Thus, we have obtained the "canonical" solution of the confluent hypergeometrical equa- 
tion in the form 



q(t) = C 1 M(a,p,t) i 



T(a) 



\ T(a-(3 + iy 



(4.16) 
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B. A canonical transformation {a, /?} —>■ 

In this case a = at and the system ( [4.3| )-( P~4"|) can be rewritten in the form 

^c-6 = d+(l--)a, t^c + b = 0, (4.17) 

a& + ^a 2 + ^"^(o' 2 - 6c) = 0. (4.18) 

We shall obtain the simplest rational solution of the inmohomogeneous undetermined system 
Q4.171 ) by taking a — 1. Then a particular solution of the obtained nonhomogeneous system 
has the form 

- p t^-\ iJr*- x l . (4.i9) 



13-/3 + 2 (3-/3 + 2 t 

From the quadratic equation Q4.18| ) we find 



(3 = (3 + 1, 7 = ^/3 -a. (4.20) 
For the generating function and canonical recurrence transformations we have 



F a ,p\ a ,(3 + x = -tV ( 2p - a qq + q 2 + (3q 2 ) , (4.21) 



d \ /— ( „ d 



y-l + j t jq=p-a~q, - \0 + t-j q = y/ '(3 - a q. (4.22) 

The recurrence relations increasing and decreasing (3 that follow from the recurrence relations 
of the handbook Q 

(3 (-1 + j) y((3) = ((3- a)y((3 + 1), - (p + tj^j y((3 + 1) = (3y((3), (4.23) 

where y{(3) = M(a,(3,t), can be obtained from Eq.( [4.22j ) by the above procedure of the 
constructing of Eq. ( |4.12j )- (|4.16| ). As a result one has 

q(t) = C 2 r(/3) M{a, (3, t). (4.24) 
T{(3 - a) 
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V. THE HYPERGEOMETRIC EQUATION 



As it was mentioned, some "radial equations" of the quantum theory are reduced to 
the hypergeometric-like equations. Therefore consider a dynamical system described by the 
hypergeometrical equation 

t(l - t)q - ((a + {3 + l)t - Qq - a(3q = 0, (5.1) 

where a, (3, ( are some parameters, as another application of the developed theory. It is the 
Lagrange-Euler equation of a dynamical system with the Lagrangian 

L = V(l - t)^-^ (> + j^/) • (5-2) 



Comparing Eq. (|ST^) with Eqs.([2.2|), (|3l3) we have: A = 0, /ii = a, /i 2 = (3, ^z — C,, 



u = Ttrrry rn = Hi-tr^\ ( 5 . 3 ) 

The system ( |5.2|) contains the three parameters a, /3, (. Therefore it is necessary to find the 
generating functions of three canonical transformations. Due to the symmetry between the 
parameters a and f3, it is sufficient to search for the generating functions of the transforma- 
tions a) {a, (3, (} — > {a, (3, (} and b) {a, f3, (} — > {a, (3, (} only. 

A. A canonical transformation {a, (3, £} — ► {a, (3, C} 



For the elementary transformation {a, (3, (} — > {a, (3, (} the system (|3.14j )- (|3.16|) can be 
rewritten the form 

(l-t r -*c-b = a-( C -- a + P-J t +1 )a, (5.4) 



(l-t)<^c + b = ^-^i3a, (5.5) 

ab = -J^f + (1 - t) a - & (bc - 7 2 ). (5.6) 

Similarly to the first case of the previous example we obtain the simplest rational solu- 
tion of the homogeneous undetermined system (|5.4j) - (|5.5|) at a = t. Then the obtained 
inhomogeneous system has the following particular solution 

In this case the quadratic equation (|5.6|) is satisfied at 
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a = a + 1, 7 = y a(a — ( + 1). 
So for the generating function we have 

i r a > /J,c|«+i,/J,C = ^ C_1 (! - [2^/a(a-C + l)gg - aq 2 - (a - ( + 1 + (3t)q 

The corresponding recurrence transformations have the form 



~2 



(5.8) 



(5.9) 



a + t— J q = y] a(a - C + 



a - C + 1 - Pt - t(l - t)— J q = ^a{a - ( + l)g. 



The standard recurrence transformation 



a + t di I V ^ = ay ( a + 1 ^ 



d 



a-C+l+Pt-t(l-t)—\q = ^Ja(a -(+ l)q, 



(5.10) 



(5.11) 



where y(a) = M(a, j3, (, t) — hypergeometric function, can be obtained by making such the 
substitution 



q{t) = A lt 



Via) 



\ r(a-c + i) 



M(a,P,(,t) 



in Eq.fCT). 



(5.12) 



B. A canonical transformation {a, (3, C} — » {a,/3,£} 
In this case the system (|3.14| )- (|3.16|) acquires the form 

^(i - «)^ c - 6 = a - (f - ° + f ~ f + 1 ) 

t c " c ~(l-t) c Hi + & = 0, 



(5.13) 



a/3 



*(1 - *) 



+ ff-C( 1 _ f )C-C( 6 _ 7 2). 



(5.14) 



We obtain the simplest solution of the system Q5.13| ) at a = 1 — t. Then a particular solution 
( |5.13| ) acquires the form 



H a — p, c 



c 



C-C-2 * t 



C-C-2 



t^- x (l-t)<- 



•C+l 



(5.15) 
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The quadratic equation ( |5.14|) leads to the following constraints 



C = C + 1, 7 = V(C-«)(0-O- 

Therefore for the generating function and the recurrence relations we find 



a,f3£\c*M+ 



(5.16) 



i = V(i - t) a+ ^ (2y/(C- ")(/?- CM + (a + /? - C)g 2 + Cg 2 ) , (5.17) 



v C-a-/3 + (l-t) 5 jg 



y/(C - a)(/3 - Qq, 



Making the substitution 



q(t) = A 2 



\ 



r(c-a)r(c-)9) 
r(C) 



M(a,/3,C,t) 



(5.18) 



(5.19) 



we reduce Eq. fl5.18|) to the relations which follow from the standard recurrence formulae 
for the hypergeometric function M(a, j3, (, t) 



(c + f|W + i) = Cv(C), 



(5.20) 



where y(C) = M(a,/?,C,0- 
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